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Abstract 

', Previous results for the prospects of B° mixing measurement in the ATLAS ex- 

Q\ ' periment at LHC are updated. The improved analysis method of the studied decay 

channels B[? — > D~7r + and Bg — > T)~af, combined with most recent values for the 
1} , branching ratios and the B° lifetime, leads to the new ATLAS sensitivity range for 

Q_i! the x s measurement: x™ ax = 42. An extensive study is done in order to estimate 

how x™ ax is influenced by the B-decay proper-time resolution of the vertex detector, 
as well as by the number of events and by the signal-to-background ratio. 

1 Measurement of x s 

CP-violation in B-meson decays can be characterized by the unitarity triangle, com- 
posed of Cabibbo-Kobayashi-Maskawa matrix elements by requiring that the matrix 
is unitary. In addition to the angles of the unitarity triangle, measurements of the 
lengths of the sides provide us with complementary information about the triangle. 
The length of the side \Vtd\ is least well-known of the sides. The mixing parameter 

is proportional to |Vf^| 2 , but inferring the value of |V^| 2 from is hampered by 
hadronic uncertainties. To large extent, these uncertainties cancel, when considering 
the ratio of the mixing parameters, x s /x& cx l^ts/V^I 2 . 

The mixing parameter x s has not been measured yet. The Standard Model pre- 
dicts it to be in the range 10-30, and the present lower limit from LEP is x s > 11.1 Q . 

Principles of the x s measurement are described in previous notes [Q, J3|. A time- 
dependent asymmetry A(t), observed in the experiment, is defined as: 

= p (++)Af - ^±4# =Dcos(x s t/r Bs ) (1) 
w dn(++)/dt + dn(+-)/dt y ' sl w 
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Figure 1: Expected amplitude of the asymmetry D = Dt a gDbackDtime( x s) as a function of 
x s for various proper time resolutions at. 



where (++) refers to B-decays in which the flavour of the B appears to be the same in 
the production and decay, (H — ) refers to cases where an oscillation of the B into its 
anti-particle is observed, and t is the proper time. D is a product of all the dilution 
factors: 
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where W is the wrong tag fraction, Ns is the number of signal events, Nb is the num- 
ber of background events^, at is the proper time resolution and r is the B^ lifetime. 
Formula |5] can be obtained using Fourier transforms 0] . 

In case of ATLAS W = 0.22 § and Nb = N s is assumed (section JHfy. The 
lifetime of the Bg is assumed to be r = 1.61 ps after Q. Figure Q shows the total 
dilution factor D estimated with the above assumptions for different values of the 
proper-time resolution (a t = 0.05, 0.07, 0.09 and 0.11 ps), as a function of x s . The 
amplitude of the cosine wave which can be observed by the experiment is reduced 

It is assumed for the formula [l] that background dNs/dt ~ exp(—t/r), with no asymmetry. 
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with increasing x s because of the finite cjj, as given by the formula || One can 
expect that for any given statistical sensitivity of the experiment, defined by the 
available signal statistics and by the background, there is a minimal value of the 
asymmetry amplitude D which still enables a reliable measurement of x s . Therefore, 
from equation ||, one expects 

x™ ax ~ — (6) 

for any objective definition of x™ ax . 



2 Definition of the ATLAS range 
2.1 Input parameters 

In order to estimate the sensitivity of the experiment for the x s measurement, we 
use as input parameters the number of signal events N$, the number of background 
events Nb, and the resolution in the proper-time of the Bg decay, ap 

Two of the Bg decay channels which can be used for the x s measurement have 
been analyzed so far: Bg — > D~7r + and Bg — > D~a^". The original work is described in 
Q and in Q, respectively. The results have been updated, for the Technical Proposal 
and on a few other occasions, because of minor changes of the cuts, updated branching 
ratios or trigger efficiencies, in particular the second level trigger sensitive to D^. It 
is therefore useful to summarize the event counting as it was done for this work. The 
calculation is explained in tables [l] and 

The number of background events for Bg — » D~7r + was also estimated in 0. 
First some potentially dangerous exclusive decays (B[J — > D~ir + , B[J — > D~7r + and 
Ab — > followed by A+ — » pK~7r + ) were checked and found not to contribute 

significantly to the background. The combinatorial background can not be precisely 
estimated because of the very large Monte Carlo statistics required. With the set 
of cuts shown in table ||, no candidate event was found in the mass window 5 < 
Mkk-kw < 6 GeV within 370,000 inclusive //X events. That gives an upper limit on 
the level of Nb <~ 0.6 • N$ within the final mass cut of ±2cr. 

A similar procedure was carried out in for the channel Bg — > Dg"a^~. The 
channels — >■ D~a^, A b — > A+7r _ followed by A+ — * pK _ 7r + 7r + 7r _ and B[J —> 
were found not to contribute significantly, as compared to the limit on combinatorial 
background, which was Nb <~ Ns- The latter was based on the Monte-Carlo sample 
of 300,000 inclusive //X events, where no event has passed cuts from table || in the 
mass window M B o - 150 MeV/c 2 < M KKn7T7T7T < M B o + 150 MeV/c 2 . 

Adding the contribution from the two channels, the total number of signal events 
is Ns = 4480 for one year of low luminosity data taking (J Ldt = 10 4 pb -1 ). A con- 
servative value of Nb = Ns was used in the following. A dedicated study of the 
impact of background on x^ ax was also done, and is described in section ||. 

The proper-time resolution at was calculated in using full GEANT simulation 
of the ATLAS tracker equipped with the strip 'B-physics' layer. The result was 
at = 0.069 ps. A slightly better value at = 0.064 ps was obtained with track level 
simulation for the second decay channel in ||. There are many other track-level 
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Table 1: Number of signal events from B° — > D s 7r + channel expected in ATLAS after 
1 year (10 7 s) of operation at 10 33 cm _2 s _1 
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B° vertex fit x 2 < 5.0 
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Track efficiency 
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results, which were produced later, and which are in a good agreement with the 
latter value. 

Replacing the strip 'B-physics' layer at R = 6 cm with a third pixel layer at 
R = 4 cm gives slightly worse performance, at = 0.068 to 0.094 ps, depending on 
the resolution assumed for pixels j/]]. Removal of the 'B-physics' layer, with the first 
tracking layer made of pixels at R = 11 cm, results in at = 0.146 to 0.154 ps. 

A value of at = 0.07 ps, well justified by the result from 0, which is the only 
one so far using full simulation, is assumed in this section. In section || the work is 
repeated for values of at ranging from 0.05 to 0.011 ps. 

A parameter of importance which was also updated for this work is the Bg life- 
time. An up to date value r = 1.61 ps used here is higher than in case of previous 
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Table 2: Number of signal events from B° — ► D s a^ channel expected in ATLAS after 
1 year (10 7 s) of operation at 10 33 cm _2 s _1 
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presentations (r = 1.54 ps, ||), adding to an increase of the predictions for x" 
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35 40 45 50 

x s 

Figure 2: An average fitted amplitude A(x s ) (left) and an average x 2 °f f ne fit (right), 
both as a function of x s , from 1000 experiments. The error bars show the RMS. 

2.2 Procedure of defining the range 

The number of signal events N$, the number of background events Nb and the 
proper-time resolution at are used as parameters in a Monte-Carlo program, which 
generates "experimental" distribution of the "measured" asymmetry A (t). The distri- 
bution then has Dback, Dtime an d Poisson fluctuations introduced with Monte-Carlo 
methods. Dilution from wrong tags Df ag is also introduced, with the wrong tag frac- 
tion W = 0.22 ||. Finally, the program generates A(t) distributions for t > 0.4 ps 
because of the analysis cut listed in tables [j] and |2[ 

Resulting 'experimental' A(t) distribution is analyzed with a method called the 
amplitude fit, recommended in A(t) distribution is fitted with the function 
Afn cos(x s t/r), where x s is a constant value and Au t is the only free parameter. The 
fit is repeated for different values of x s , giving an Afn{x s ) distribution. 

The amplitude fit method can be described as a variant of the Fourier transforma- 
tion, a 'cosine transformation'. It has the power of Fourier analysis for the periodical 
signals, but one introduces known information about the phase in addition - we are 
looking explicitly for a cosine wave in the A(t) distribution. 

Figure || shows the properties of the transformation. An average fitted amplitude 
is shown as a function of x s after 1000 experiments which had the true x s = 42. An 
average x 2 value of the fit is also shown as a function of x s and the error bars show 
the RMS calculated from 1000 'experiments'. It can be observed, that Afn(x s ) is 
more useful to find the true value of the oscillation frequency. 

The value of x s giving the highest Aja is the x s measurement of the experiment. 
The peak in the fitted amplitude has its natural width, which can be seen in Fig. [2| 
(more explanation can be found in Q). An experiment is called 'successful' here 
when the measured x s value is within that width from the 'true' value defined in the 
Monte-Carlo. 
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Figure 3: Probability of a successful measurement, estimated using 1000 "experiments" 
with nominal ATLAS parameters: a t = 0.07 ps and / Ldt = 10 4 pb _1 . 



With increasing 'true x s f the amplitude of the oscillation seen in the A(t) distri- 
bution decreases, as shown in figure [l|. When a limit is passed the amplitude fit no 
longer enables to distinguish the right oscillation frequency form the noise generated 
by the statistical fluctuations in bins. 

The probability of the experiment success calculated from 1000 simulated 'exper- 
iments' is shown in figure | as a function of x s . The limit x™ ax is defined as the 
highest value of x s for which such probability of an experimental success is above 
95%. 

For ATLAS, with the parameters listed in section 2A the limit turns out to be 
x max _ 42. A time-dependent asymmetry distribution, its transformation with the 
amplitude fit procedure and the results of 1000 simulated experiments are shown in 
Fig.|. 

For the correct frequency an average amplitude should be equal to the product 
of all dilution factors D (see formula [l]) . Figure ^ shows the distribution of the fitted 
amplitudes. An average value (53.15 ±0.49) • 10 -3 is in agreement with the prediction 
from formulas || to | equal D = 52.85 • 10 -3 . 
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cr t = 0.07 ps, 10 4 pb" 1 
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Figure 4: x s = 42, "nominal" parameters: a t = 0.07 ps and / Ldt = 10 4 pb -1 ; (a) time- 
dependent asymmetry distribution (single experiment), (b) amplitude A(x 8 ) fitted to the 
distribution shown in (a), (c) distribution of x s giving the highest value of A(x s ) when the 
experiment was repeated 1000 times. 
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Figure 5: Distribution of the amplitudes fitted for x s = 42 for 1000 'experiments' (generated 
with the 'true' x s = 42). Dark histogram shows amplitudes for 31 'experiments' which 'fail'. 



3 Dependence of the range on at and on inte- 
grated luminosity 

In order to evaluate the importance of the proper time resolution at and of the inte- 
grated luminosity, the analysis was repeated for different values of these parameters. 
All combinations given by at values of 0.05, 0.07, 0.09 and 0.11 ps, together with 
J Ldt values of 0.5, 1.0, 2.0, 5.0 and 10.0 x 10 4 pb" 1 were tried. 

Deteriorating the proper time resolution form the 'nominal' 0.07 ps to 0.11 ps 
reduces the x™ ax to 27, as could be expected from formula ||. The corresponding 
summary plot is shown in figure ||. 

If at has the nominal value and the integrated luminosity is reduced to 0.5 • 
10 4 pb -1 the x™ ax is less reduced: x'f ax = 38 (see figure 0). 

A full scan of the two parameters also gives a set of values: at = 0.05 ps and 
/ Ldt = 10 5 pb _1 . Such a combination gives x™ ax = 77 (Fig. §). 

A summary plot showing the x™ ax versus at for various integrated luminosities 
is shown in figure []. The expected dependence given by the formula [], shown with 
solid lines in the figure, is well confirmed. 

The same summary in another projection, x™ aa: versus / Ldt for different values 



of at, is shown in figure 10. It is found, that the ATLAS range for the x s measurement 
can be parameterized in the following way: 
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where at is given in ps and / Ldt in units of 10 4 pb 1 
show the prediction of the above formula. 



The solid lines in figure IK 
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cj t = 0.1 1 ps, 10 4 pb 
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Figure 6: x s = 27, o t = 0.11 ps and / Ldt = 10 4 ptr 1 ; (a) time-dependent asymmetry 
distribution (single experiment), (b) amplitude A(x 8 ) fitted to the distribution shown in 
(a), (c) distribution of x s giving the highest value of A(x s ) when the experiment was 
repeated 1000 times. 
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cr t = 0.07 ps, 0.5*1 4 pb" 1 
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Figure 7: x s = 38, a t = 0.07 ps and / Ldt = 0.5 • 10 4 pb -1 ; (a) time-dependent asymmetry 
distribution (single experiment), (b) amplitude A(x s ) fitted to the distribution shown in 
(a), (c) distribution of x s giving the highest value of A(x s ) when the experiment was 
repeated 1000 times. 
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Figure 8: x s = 77, a t = 0.05 ps and J Ldt = 10 • 10 4 pb" 1 ; (a) time-dependent asymmetry 
distribution (single experiment), (b) amplitude A(x s ) fitted to the distribution shown in 
(a), (c) distribution of x s giving the highest value of A(x s ) when the experiment was 
repeated 1000 times. 
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Figure 9: Sensitivity range of the ATLAS experiment for the x s measurement, as a function 
of the proper-time resolution a t , for various integrated luminosities. The points come from 
Monte-Carlo calculations, the solid lines show l/a t dependence. 



4 Dependence on signal-to-background ratio 

The scan of parameters at and / Ldt described in the previous section was done for 
the background-to-signal ratio of 1. By increasing the number of background events 
such that Ns/(Nb + Ng) becomes equal 0.1 we reduce x™ ax to 31 for the 'nominal' 
parameters at = 0.07 ps and / Ldt = 10 4 pb _1 (see figure 11). 

A summary plot showing the dependence of x™ ax on background is shown in figure 
|T^ . Ns was kept constant at its predicted value for / Ldt = 10 4 pb _1 (N$ = 4480, 
see section ^J), while the number of background event was varied between Nb = 
and Nb = 9 x Ns- It can be observed, that the dependence of x™ ax on background 
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Figure 10: Sensitivity range of the ATLAS experiment for the x s measurement, as a func- 
tion of the integrated luminosity, for various proper-time resolutions a t . The points come 
from Monte-Carlo calculations, the solid lines show the x™ ax ~ (2.90/<r t )- (1+0.15 ln(/ Ldt)) 
parametrisation, where a t is in ps and / Ldt is in units of 10 4 pb^ 1 . 



is not strong. It should be noted however, that background had no asymmetry in the 
simulation and a simple exponential dependence dNs/dt ~ exp(— t/r) was assumed. 
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cr t = 0.07 ps, 10 4 pb _1 , N S /N T = 0.1 
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Figure 11: x s = 31, a t = 0.07 ps and / Ldt = 10 4 pb _1 , N B /N S = 1.0; (a) time-dependent 
asymmetry distribution (single experiment), (b) amplitude A(x 8 ) fitted to the distribution 
shown in (a), (c) distribution of x s giving the highest value of A(x s ) when the experiment 
was repeated 1000 times. 
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Figure 12: Sensitivity range of the ATLAS experiment for the x s measurement, as a func- 
tion of the signal content of the sample, for a t = 0.07 ps and / Ldt = 10 4 ptr 1 . 
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5 Conclusions 



It is justified theoretically, after that the sensitivity range x™ ax is proportional to 
l/cr t (where at is the resolution in the B-decay proper time), for any given integrated 
luminosity and signal-to-background ratio, independently of the method of defining 
the sensitivity range. 

With the methods described in section |2| it is found that x™ ax = 42 for nominal 
ATLAS parameters. The Monte-Carlo study done for various integrated luminosities 
and for different at values shows, that ATLAS sensitivity range can be parameterized 
as x™ ax ~ (2.90/oi) • (1 + 0.15 ln(J Ldt)), where a t is in ps and / Ldt is in units of 
10 4 pb _1 . It is also found, that the dependence of x™ ax on background is relatively 
weak. 
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